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, We analyze a class of multicomponent nonlinear Schrodinger equations (MNLS) related to the symmetric BD.I- 
type symmetric spaces and their reductions. We briefly outline the direct and the inverse scattering method for 
^ , the relevant Lax operators and the soliton solutions. We use the Zakharov-Shabat dressing method to obtain the 
' two-soliton solution and analyze the soliton interactions of the MNLS equations and some of their reductions. 
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1. INTRODUCTION 



. Bose-Einstein condensate (BEC) of alkali atoms in the F — 1 hyperfine state, elongated in x direction and 
confined in the transverse directions y, z by purely optical means are described by a 3-component normalized 
^ . spinor wave vector = ($1, $0, t) satisfying the nonlinear Schrodinger (MNLS) equatiorP see 

alsoPHSl 

> ■ idt^x + 52*1 + 2(|$i|2 -t- 2|$on$i + 2$*_i$^ = 0, 

i9t$o + a>o + 2(|$_i|2 + |$o|' + |$i|')*o + 2$S*i$-i -0, (1) 
idt^-x + + 2(|$_i|2 + 2|$on$-i + 2$t$g = 0. 



spinor BEC with F = 2 for rather specific choices of the scattering lengths in dimensionless coordinates takes 
the form.*^ 



idt<^±2 + 9xx$±2 + 2(*, **)$±2 - (2$2$-2 - 2$i$_i + $g)$^2, 

idt^±i + a^.^ii + 2(*, #*)$±i - (2$2$-2 - 2$i$_i + (2) 

2\ 



[X] , idt^o + + 2(*, **)$o - (2$2$-2 - 2$i$_i + 



Both models have natural Lie algebraic interpretation and are related to the symmetric spaces BD.I ~ 
SO(n -I- 2)/S0(n) x S0(2) with n = 3 and n — 5 respectively. They are integrable by means of inverse scattering 
transform method.l^''^'i2l Using a modification of the Zakharov-Shabat 'dressing method' we describe the soliton 
solutiontP^i^ and the effects of the reductions on them. 



Sections 2 and 3 contain the basic details on the direct and inverse scattering problems for the Lax operator. 
Section 4 outlines the effects of the algebraic reductions of the MNLS. In Section 5 using the Zakharov-Shabat 
dressing method we derive the one- and two-soliton solutions of the MNLS and discuss their properties. Section 
6 is dedicated to the analysis of the soliton interactions of the MNLS. To this end we evaluate the limits of the 
generic two-soliton solution for t — ^ ±00. As a result we establish that the effect of the interactions on the soliton 
parameters is analogous to the one for the scalar NLS equation and consists in shifts of the 'center of mass' and 
shift in the phase. 
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2. THE METHOD FOR SOLVING MNLS WITH F = 1 AND F = 2 

MNLS equations for the BD.I. series of symmetric spaces (algebras of the type so{n + 2) and J dual to ei) have 
the Lax representation [L, M] = as follow^E"i2l 



(3) 



L-)p{x, t, A) = i^j^V + U{x, t, X)ip{x, t, A) = 0, Mip{x, t, A) = idtip + V(x, t, X)i:{x, t, A) = 0, 

U{x,t,X) ^Q{x,t)-XJ, V{x,t,X) ^ Voix,t) + XVi{x,t) ~ X^J, 

Vi{x,t) = Qix,t), Vo{x,t) = ^adji^ + ^ [adj^Q,Q{x,t)] . 

where ad jX — [J, X] and ad is well defined on the image of ad j in g; 

/ ,f \ 

Q=\P* soq \ , J = diag(l,0,...0,-1). (4) 
V pUo J 

The vector q for F ^ 1 (resp. _F — 2) is 3- (resp. 5-) component and has the form 

g = ($i,$o,*-i)^, q = ($2,*i,*o,$-i,*-2)^, (5) 
and the corresponding matrices sq enter in the definition of so{2r + 1) with r = 2 and r = 3: 

2r+l / 1 \ 

Xeso(2r + l), X + ^o^^^o = 0, 5o = = -sq . (6) 

s=i \ 1 / 

By i?sp above we mean 2r + 1 x 2r + 1 matrix with matrix elements {Esp)ij ~ SsiSpj. With the definition of 
orthogonality used in (|6]) the Cartan generators Hk = Ek.k^ En+3-k,n+3-k are represented by diagonal matrices. 

If we make use of the typical reduction Q = (or p* — q) the generic MNLS type equations related to 
BD.I. acquire the form: 

iqt + qxx + 2{q\q)q- {q, SQq)Soq* = 0. (7) 
The Hamiltonians for the MNLS equations ([7]) are given by 



dx (^{d^q-\d.,q) - {q\q)^ + i|(g^,,o9)p) 



ffMNLS - / dx[ {d,q\d.,q) - {q\q)^ + -\{f, s^q)? ■ (8) 



3. THE DIRECT AND THE INVERSE SCATTERING PROBLEM 
3.1 The fundamental analytic solution 

We remind some basic features of the inverse scattering theory for the Lax operators L, see.'^'i^! There we 
have made use of the general theory developed irP^HlSlIHl g^j^^j ^jjg references therein. The Jost solutions of L are 
defined by: 

lim 0(a;, i, A)e^^'^"^ = 1, lim i/;(x, t, A)e'^"'^ = 1 (9) 



x—^ — oo 



and the scattering matrix T{X,t) = ip ^4>{x,t,X). The special choice of J and the fact that the Jost solutions 
and the scattering matrix take values in the group SO{2r + 1) we can use the following block-matrix structure 
of r(A,i) 



r(A,i) 



/ 771+ -b-^ cr 

b+ T22 -soB- I , T{X,t) = 



I 





B-^ 




-B+ 


T22 


sob~ 






777,^ 



(10) 



where 6^(A,i) and B^{X,t) are 2r — l-component vectors, T22(A) is 2r — 1 x 2r — 1 block matrix, and mf{\), 
and cf{X) are scalar functions. Such parametrization is compatible with the generalized Gauss decompositions 
of T(A) which read as follows: 



T(A,t) = TJD+Sj, T{X,t) = T+DjSj, 

TJ = e±('^*'^f ), = e±(^~*'^?), = diag ((mf )±i,m±, (m±)Ti) 



(11) 



where 



fc=i 



(12) 



fc=l 



and similar expressions for (^T^y '^^) ' '^^^ functions mf and n x n matrix- valued functions are analytic for 
A e C±. We have introduced also the notations: 

TO^ ™1 '^1 ^7^1 

There are some additional relations which ensure that both r(A) and its inverse T(A) belong to the orthogonal 
group SO{2r + 1) and that T{X)f{X) = 1. 

Next we introduce the fundamental analytic solution (FAS) x*(a;,t, A) using the generalized Gauss decom- 
position of r(A,i), seePIiSMl 

X^{x, t, A) = ^{x, i, X)Sf{t, A) = ib{x, t, X)Tf{t, X)Df{X), (13) 

This construction ensures that ^^(a;,A) — x^(a;, A)e*'^'^^ are analytic functions of A for A e C±. If Qjx, t) is a 
solution of the MNLS eq. ([7]) then the matrix elements of r(A) satisfy the linear evolution equationai^] 

^^±A26±(t,A) =0, i^±X^B^{t,X)=Q, 
at at (14) 

dt dt 

Thus the block-diagonal matrices D^{X) can be considered as generating functionals of the integrals of motion. 
The fact that all (2r — 1)^ matrix elements of mJ(A) for A € C± generate integrals of motion reflect the 
superintegrability of the model and are due to the degeneracy of the dispersion law of ([7]). Note that i'j (A) 
allow analytic extension for A G C± and that their zeroes and poles determine the discrete eigenvalues^ of L. 

3.2 The Riemann-Hilbert Problem 

The FAS for real A are linearly relatecP^ 

X+{x,t,X) = x-{x,t,X)G^,j{X,t), Go,j{X,t) = Sj{X,t)Sj{X,t) (15) 
Eq. (fT5|) can be rewriten in an equivalent form for the FAS £,^{x, t, A) = X^{xi t, X)e^^'^^: 

t^+Qix)^^{x,X)-X[J,^^{x,X)]=0, (16) 
dx 



and the relation 



lim ^±(x,t, A) = 1, (17) 

A— >oo 



Then these FAS satisfy the RHP's 



(18) 



Obviously the sewing function Gj{x,X,t) is uniquely determined by the Gauss factors Sf{X,t). In addition 



Zakharov-Shabat's theorenP^'^l g^;ates that Gj{x, X,t) depends on x and t in the way prescribed above then the 
corresponding FAS satisfy the hnear systems (fTB]) . 

If we have solved the RHP's and know the FAS £,^{x,t,X) then the formula 

Q{x,t)^ lim X(j -^+{x,t,X)ji+{x,t,X)] , (19) 
allows us to recover the corresponding potential of L. 



4. REDUCTIONS OF MNLS 

Along with the typical reduction Q — mentioned above one can impose additional reductions using the 
reduction group proposed by Mikhailov.123 They are automatically compatible with the Lax representation of 
the corresponding MNLS eq. Below we make use of two Z2-rcductionsP2l 

1) CiU''{x,t,X*)G^^ = U{x,t,X), GiV''{x,t,X*)G{^ =V{x,t,X), 

2) C2U^{x,t,X)G2^ ^ -U{x,t,X), C2V'^{x,t,X)G2^ = -V{x,t,X), 

where Ci and C2 are involutions of the Lie algebra so{2r + 1, i.e. Cf = 1. They can be chosen to be either 
diagonal (i.e., elements of the Cartan subgroup of S0{2r + 1)) or elements of the Weyl group. 

The typical reductions of the MNLS eqs. is a class 1) reduction obtained by specifying Ci to be the identity 
automorphism of g; below we list several choices for Ci leading to inequivalent reductions: 



(20) 



where 

and e," 



la) 


Ci = a, 


p{x) 




lb) Ci = 


Ki, 


p{x) 


= Koiq*{x), 




Ic) 


Ci = Se^ , 


p{x) 


Ko2q*{x), 


Id) Ci = 




p{x) 




(21 


2e) 


C2 = K4, 


q{x) 


= ~Ko4:Saq{x), 


p{x) = 


-Ko4Sop{x), 










— block-diag 




, 1), Koi = 


diag(ei, . . . , 


-li 1, ^r-1, ■ ■ 


,ei), 


J = 1,2,3, 


(22 



±1. The matrices Kq2, Kq^ and are not diagonal and may take the form: 



K02 = 



K02 = 


































-1 


i) 








( 













1 















1 























-1 \ 




( 



















1 
























1 











-1 








Ko3 - 










-1 

















1 












1 











-1 











J 




[ 


-1 











/ 



(23) 



V 

Each of the above reductions impose constraints on the FAS, on the scattering matrix T{X) and on its Gauss 
factors Sf{X), Tf{X) and Dj (A). These have the form: 

(5+(A*))t = Kr^S-{X)K, (r+(A*))t = KJ^f-{X)K, (D+iX*))^ = KJ^b-{X)K, 

where the matrices Kj are specific for each choice of the automorphisms Ci, see eqs. (HU, 1211) ■ 
In particular, from the last line of and we get: 

(m+(A*))*=™r(A), 
+ 



(24) 



and consequently, if toJ'(A) has zeroes at the points A^, then (A) has zeroes at: 



(25) 
(26) 



5. SOLITON SOLUTIONS 

Let us now make use of one of the versions of the dressing methocP^E^ -vvhich allows one to construct singular 
solutions of the RHP. In order to obtain 7V-soliton solutions one has to apply dressing procedure with a 2A^-poles 
dressing factor of the form 



u[x 



The 7V-soliton solution itself can be generated via the following formula 

N 

Qjv,s(a:^) = ^[J,A,(x)+Bfc(a;)]. (28) 

fc=i 

The dressing factor u(x, A) must satisfy the equation 

idxU + Qn^sU — X[J,u] ~ (29) 

and the normalization condition lim^^oo w(x, A) = 1. The construction of u(x,X) £ SO{n + 2) is based on an 
appropriate anzatz specifying the form of its A-dependenc^^HUll 

The residues of u admit the following decomposition 

Ak{x) = Xk{x)Fl{x), Bkix) - Ykix)Gl{x). 

where all matrices involved are supposed to be rectangular and of maximal rank s. By comparing the coefhcients 
before the same powers of A - A^ in (EH) we convince ourselves that the factors Fk and Gk can be expressed by 
the fundamental analytic solutions Xoi^^ ^) follows 

F^i^) = Flo[xU^^^t)r\ Gl{x) = Gl,[x^{x,\l)]-\ 
The constant rectangular matrices i^l-.o and Gk.o obey the algebraic relations 

Fk,o^oFk,o — 0, G^qS'oG/co — 0. 
The other two types of factors and are solutions to the algebraic system 



bol^k - AfeOfc + 2^ — + 2^ — -p, 

T T (30) 

boGk = 2^ — p — + Yklik + 2^ — 

I l^k ^ ^k 



The square s x s matrices ak{x) and Pk{x) introduced above depend on Xo a-^d Xo ^iid their derivatives by A 
as follows 

ak{x) = -F^j^[xt{x,Xl)]-^dxxt{x^K)SoFo.k+aQ^k, ^^^^ 
lik{x) = -G'^f.lxQ [x, dxxn {x,X^)SQGo^k + Po,k- 

Below for simplicity we will choose F^. and Gk to be 2r + 1-component vectors. Then one can show that 
ctk — Pk = which simplifies the system ([5(1)) . We also introduce the following more convenient parametrization 
for Fk and Gk, namely (see eq. (|55|) ): 

Fkix,t) ^ So\nk{x,t)) = \ -72so% , Gk{x,t) = \nl{x,t)) ^ i V^uok* I, (32) 



where uok are constant 2r — l-component polarization vectors and 

Zj =iyj{x + 2fijt) + Coo, <i>j = IJ-jX + in] - v'^j)t + 5qq, 
{nj{x,t)\So\nj{x,t)) =Q, or {vqjSqVqj) = I. 

With this notations the polarization vectors automatically satisfy {nj{x,t)\SQ\nj{x,t)) = 0. 
Thus for A'' = 1 we get the system: 

(A+-Ar)|ni) (A+ - K)So\nt) 

which is easily solved. As a result for the one-soliton solution we get: 

gis = - ^ ' (e-^^5o|%) + e^M^i», Ai = cosh(2zi) + (4|%). (35) 

For n = 3 we put vok = |i'o/c|e"°'° get: 
*is;±i = - ^/^l'^°^'^'^°l^3|(A^-Al) ^_,^^^,^^3 ^^^^^^^^ ^ ^^^^ ^^^^^^^^ - ^sinh(2l T Coi) sin(ai3)) , 
^^^^^ ^ _ V^koi;2|(A+-Ar) ^_,^, ^^.^j^^^ ^^^^^^^^ + icosh^i sin(ao2)) , (36) 

As = ^(ao3 - aoi), Coi = \ ' "13 = ^("03 + aoi), 

Note that the 'center of mass' of <l>is;i (rcsp. of $is;_i) is shifted with respect to the one of $is;o by Coi to the 
right (resp to the left); besides |$is;i| — \^is;-i\, i-c. they have the same amplitudes. 

For n = 5 we put i^ok = |t'ofc|e"'"= and get analogously: 

^^^^^^ ^ _ V2Woi:i'^oi,,\{Xt ^K)e-i<^M. (eosh(^i t Coi)cos(ai5) - isinh(^i t Coi) sin(ai5)) , 
^1 

*is;±i = V^I'^oi-2'^°i^^l(^t ^Ale-»^i±'fe4 (cosh(zi T C02) cos(a24) - »sinh(zi T Coi) sin(a24)) , 
^1 

^^^^^ ^ _ ^l'^oi;3|(Ai^-Ar) ^_,^, ^^^^^^^ ^^^^^^^^ _ .^^^^^ ^.^^^^^^^ ^ (37^ 

As = ^(ao5 - aoi), Coi = \ ^\iAn^i\ ' "is = ^("05 + "oi), 

/324 = ^(Q!04 - "02), C02 = :^ In 1^^^^ , "24 = ^(ao4 + 002), 

Similarly the 'center of mass' of <i?is;2 and $is;i (resp. of <&is;-2 and $is;-i) are shifted with respect to the one 
of $is;o by Coi and C02 to the right (resp to the left); besides |$is;2| = |^is;-2| and |$is;i| = |^is;-i|- 
For N = 2 we get: 

\ni{x,t)) 

\n2{x,t)) 
So\nUx,t)) 
So\n2{x,t)) 



X2{x,t)f21 


Fi(.t,<)kii 


1 l2(a;,t)K2i 


X2 - xt 


Ar-A+ 


A2 - A^ 


Xl{x,t)h2 


Fl(x,t)Kl2 


1 Y2{x,t)K22 


Xi - X2 


X^ - X2 


A2 — 


Xi{x,t)Kxi 


1 X2{x,t)Kii 

xt - Xi 


1 Y2{x,t)f^, 

X2 ^ '^1 


xt " xt 


Xi{x,t)K2l 


1 X2{x,t)n22 


Yiix,t)f^2 


xt - A2 


xt ~ X2 


Xi - X2 



(38) 



where 



fkj{x,t) = e'"-'^-'^''"'-'*'^'' + 5^^-^'=+*^'=-'^^) - 2 (j7(ffeSo!7oj) , 



In other words: 



MX 



( 

X+ \ + 
K21 



Aj -Aj 



K22 



A--A+ 

K12 
Aj -A2 



/,* 



\ Aj — Aj Aj — Aj Aj — Aj 

We can rewrite M in block-matrix form: 

Mu M12 



A2 -A]^ 

«22 
^2 ^-''2 

A--A- 










( \ni) \ 








\n2) 












\ ^2 j 







The inverse of M is given by 



-22 



M90 = M 



11; 



M21 = -Mfa, 



12 



A+-A+ V -1 



1 



K21 



M12 = 



Aj A2 — A"!*" 

K12 K22 



(Mil -Mi2M^iM2i)-i 
-(M^i - M2iMiiMi2)-iM2iMii 

One can check by direct calculation that: 

Mil - Mi2MtiM2i = 



-(Mil - Mi2MtiM2i)-^Mi2M^ 
(Mti -M2iMiiMi2)-i 



/ 



12 



Mil -M21M11M12 
Z 



I/12P 



A2 - 

/12 

AJ - Aj^ 

K12K21 



1 

-1 

1 

-1 

K11K22 



IA2+-A+I2 |A+-ArP Mz.2 



Finally we get: 



M" 



1 



/ 



V 





/1* 



Aj -A2 
K22 

_ K12 

a+-a: 



fl2 

\--\- 


_ K21 

Aj^-Aj 



A+-A,- 



/12 

A+-A+ 



K12 \ 

A2 -A;^ 
_ Kll 

112 
^1 ^^2 





From eqs. (|40|) and (|44|) we obtain: 



1^2) 

|i^2> 



1 

^ Va 
1 



/r2 



1 



1 ^ ^^2 
^12 
^1 ^ ^2 
K22 



'*22 I *\ 

^2 ^ -^2 



'«12 



\ni) 



K21 



Z \ A^ — A9 



-m - 



1 
z 



K12 



AJ - Aj 



"1 



Ai - A2 
K21 
- A2 
^11 

A+-Ar 



■So\nl) 



A2 - Aj 



-5'o|"2> 



A+-AI 
/12 



-So\nl 



5, 



A+-A+ 
AJ - Ai 



■0|'^2, 



-So\nl 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



Inserting this result into eq. (pS)) we obtain the following expression for the 2-soliton solution of the MNLS: 



Q2s{x, t) = [J, ^1 + Si + A2 + B2] = ^ [J, t) - SoC^ix, t)Sol 



C{x,t) = ,+^^\- 1^1) K I - ,+^^\- l'^l)(4l - T+~^^l'^2>(?4l + ,+'^",- l"2)(»2l (-46) 



—^^\ni){n2\So - ■T+-^^So\n*2){n\\. 
^1 ^ ^2 -^1 ^ ^2 



At the end of this section we note that the effect of the reductions (|20|) ~([2T 1) consists in constraining the 
polarization vectors. For the reduction 2e) we get 

vok = KoiVok (47) 

In particular, for n = 3 and for K04 = — 1 we have qi = qs, and 52 arbitrary. This reduction of eq. ([T]) is also 
important for the BEC- From (|T7|) we find i/qi — 1^03. The effect of this constraint is that for the one-soliton 
solution we get $is;i = $is;3- 

Our next remark followinJ^H is that this reduction applied to the F = 1 MNLS ([!]) leads to a 2-component 
MNLS which after the change of variables 

qi^^{wi+iw2), q2 = ^{wi - iw2), (48) 

leads to two disjoint NLS equations for wi and W2 respectively. 

It is only logical that applying the constraint vqi = ^'os the explicit expression for the one-soliton solution 
(1551) simplifies and reduces to the standard soliton solutions of the scalar NLS. 

6. TWO SOLITON INTERACTIONS 

In this section we generalize the classical results of Zakharov and Shabat about soliton interaction^^!] to the class 
of MNLS equations related to BD.I symmetric spaces. For detailed exposition see the monographs .'^^"^^ These 
results were generalized for the vector nonlinear Schrodinger equation by Manakovf^see also.^^^^^ The Zakharov 
Shabat approach consisted in calculating the asymptotics of generic iV-soliton solution of NLS for t — > ±00 and 
establishing the pure elastic character of the generic soliton interactions. By generic here we mean iV-soliton 
solution whose parameters = /ife ± ivk are such that /i^ ^ /ij for k ^ j. The pure elastic character of the 
soliton interactions is demonstrated by the fact that for t ±00 the generic A''-soliton solution splits into sum 
of N one soliton solutions each preserving its amplitude 2vk and velocity //fe. The only effect of the interaction 
consists in shifting the center of mass and the initial phase of the solitons. These shifts can be expressed in terms 
of only; for detailed exposition see.^i^ 

We start with the simplest non-trivial case. Namely we use the 2-soliton solution derived above and calculate 
its asymptotics along the trajectory of the first soliton. To this end we keep zi{x, t) fixed and let r = 2:2 — zi tend 
to ±00. Therefore it will be enough to insert the asymptotic values of the matrix elements of M for r ±00 
and keep only the leading terms. That gives: 

K22 = 



K12 



K21 



e2^exp(i^2Zi/j^i) + 261, 


for r cx), 






exp(-z^2Zi/i^i) + 261, 


for T —00, 






exp((l + V2/vi)zi + i(0i - (, 


^'2)) + 0(l), 


for 


r — cx). 


e^^ exp(-(l -1- i'2/i'i)zi - i{4>i 


-02)) +0(1), 


for 


r — —00 


exp((l + !^2/'^i)^i - *(0i - <; 


^'2)) + 0(l), 


for 


r 00, 


e"^ exp(-(l -1- V2lvi)zx + i{(i>\ 


-02)) +0(1), 


for 


T — > —00 


exp(-(l - 1^2/1^1)2:1 + i{4>\ - 


-02)) + 0(1), 


for 


T 00, 


e~''exp((l - V2lvx)zx - i((/)i - 


02)) + 0(1), 


for 


r — )• —00, 



(49) 



After somewhat lengthy calculations we get: 



lim q2s{x,t) 

T—>00 

lim q2s{x,t) 



iV2z/ie-'("^i-«+) (e-^i-''+so|i7oi) +e^'+''+Ki)) 
cosh(2(^;i +r+)) + {u^^, uqi) 
i\/2z/ie-*('^i+«+) (e-^i+''+So|i7oi) + e^i-''+|!^i)) 
cosh(2(2;i - r+)) + (i^iJi, i^oi) 



(50) 



where 



r+ = In 



4 



Xf - X2 



a+ = arg 



A+ - A2 



In other words the 2-soliton interaction for the MNLS eqs. related to the BD.I symmetric spaces is the same 
as the one of the scalar NLS. Again we have that for large times the 2-soliton solution splits into sum of 1-soliton 
solutions with shifted center of masses and phases and the value of these shifts r+ and a+ are independent on 
the number of components of MNLS. It will be interesting to check whether the AT-soliton interactions consist 
of sequence of elementary 2-soliton interactions and the shifts are additive. 



7. CONCLUSIONS AND DISCUSSION 

Using the Zakharov-Shabat dressing method we have obtained the two-soliton solution and have used it to 
analyze the soliton interactions of the MNLS equation. The conclusion is that after the interactions the solitons 
recover their polarization vectors t'ofe, velocities and frequency velocities. The effect of the interaction is, like in 
for the scalar NLS equation, shift of the center of mass zi — >■ zi + r+ and shift of the phase ^ (pi + a+ . Both 
shifts are expressed through the related eigenvalues A^ only. 

The next stop would be to analyze multi-soliton interactions. Our hypothcslK is that each soliton will acquire 
a total shift of the center of mass that is sum of all elementary shifts from each two soliton interactions. Similar 
result is expected for the total phase shift of the soliton. Proofs of these facts will be published elsewhere. 
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